Probing biological light-harvesting phenomena by optical cavities 
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We propose a driven optical cavity quantum electrodynamics (QED) set up aimed at directly probing en- 
ergy transport dynamics in photosynthetic biomolecules. We show that detailed information concerning energy 
transfer paths and delocalization of exciton states can be inferred (and exciton energies estimated) from the 
statistical properties of the emitted photons. This approach provides us with a novel spectroscopic tool for the 
interrogation of biological systems in terms of quantum optical phenomena which have been usually studied for 
atomic or solid-state systems, e.g. trapped atoms and semiconductor quantum dots. 



Introduction- Following a series of fascinating non- 
linear optical experiments with light-harvesting complexes 
(LHCs) involved in natural photosynthesis |1-5|, the anal- 
ysis of exciton dynamics in biological systems has attracted 
a rapidly growing interest from different scientific communi- 
ties. It has been demonstrated that the essentially lossless ex- 
citation energy transfer (EET) in these structures does not fol- 
low the weak-coupling Forster theory 1 6 1 where the excitation 
energy incoherently hops between single molecules. Indeed, 
exciton states in LHCs are frequently delocalized over several 
molecules. Several theoretical studies have analyzed EET in 
such pigment-protein complexes and demonstrated that an in- 
tricate interplay between environmental noise (dephasing) and 
quantum coherence is crucial to explain the remarkably high 
EET efficiency (above 95%) from the peripheral antenna, col- 
lecting external light, to a reaction center, where electronic 
energy is transformed into chemical form |7 - 17 |. 

On the other hand, the ability to control light-matter interac- 
tion for single atoms ifTSl or quantum dots 119 1 in high-quality 
optical cavities, e.g. photonic crystals, micropillars, and mi- 
crodisk resonators, has recently allowed the theoretical and 
experimental analysis of fundamental quantum phenomena 
arising at a single photon level |20|. In particular, quantum 
optical phenomena such as the vacuum Rabi splitting, single- 
atom lasing, photon squeezing, anti-bunching, and matter- 
light entanglement have been explored in the strong-coupling 
limit, i.e. when the exchange energy rate between the cavity 
mode and the trapped system is larger than the involved de- 
cay processes ||2TJ|22l. In this regime, the cavity and trapped 
system can no longer be described as separate subsystems, 
i.e. there is a splitting between energy eigenstates into differ- 
ent manifolds (dressed states or polaritons associated to the 
Jaynes-Cummings model). Polariton physics has had a recent 
surge in interest in the context of coupled cavity arrays 123] - 
|25]| and in Bose-Einstein condensation of exciton polaritons 
(2611271. A different, yet related, lasing process 1281 was also 
observed in organic molecular crystals in cavities 1291 . 

Here, we present the extension of the rich physics of cavity 
quantum electrodynamics to the investigation of the internal 
excitation energy transfer dynamics of bio-molecules. In par- 
ticular, we propose theoretically an experimental scheme in 




FIG. 1: Artistic impression of the proposed setup: a sample of light- 
harvesting complexes is confined in an optical cavity (resonance fre- 
quency cjc) with a quality factor Q. The sample is subjected to an 
orthogonal driving laser field of frequency uji . Photons leaking out 
of the cavity are collected by a photon detector. 



which a sample of LHCs, driven by an external laser field, is 
trapped inside an optical cavity. Proof-of-principle demon- 
strations for solution-phase cavity-enhanced spectroscopy 
were reported in |30|, where absorption measurements with 
Bacteriochlorophyll a (BChla) in solution were performed. In 
this Letter, the main goal is to detect quantum phenomena 
associated with the emitted quantized light |31|, in order to 
probe biological coherence in photosynthetic complexes, in a 
different and complementary way with respect to 2D nonlin- 
ear electronic spectroscopy used so far in Refs. [iTl-[4ll36l. 

We show that the LHC exciton structure and energy transfer 
paths can be mapped onto the statistical properties of cavity 
photons such as mean photon number and second-order corre- 
lation functions and later measured through a light mode leak- 
ing out of the cavity. Furthermore, we quantify the amount of 
quantum correlations between the cavity mode and the biolog- 
ical system and discuss the possibility of polariton formation 
for a single LHC in a cavity. 

Model- In order to probe the exciton structure of biolog- 
ical systems, we essentially propose a pump-probe scheme in 
which the probe field is substituted by the cavity mode. The 
sample (system) is assumed to be confined inside an optical 
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cavity with resonance frequency Uc, while energy is injected 
into the system via an external laser field (pump) with fre- 
quency uji . While the pump field can be pulsed or continuous 
wave (CW), the probe field (cavity mode) is always CW. We 
will consider the case where the pump field is turned on at a 
certain time, which allows us to explore both the stationary 
response to the CW field as well as the transient excitation 
dynamics. Moreover, for the sake of simplicity, we consider 
a collinear polarization of the laser field and the cavity mode. 
The extension of the model to a pulsed excitation and cross- 
field polarizations is straightforward. The LHC exciton dy- 
namics is characterized by monitoring the mean photon num- 
ber and the second order coherence function ^^^^(0) of the 
light leaking out of the cavity (probe). The coherent interac- 
tion between the cavity mode and the LHC chromophores can 
be written as 

He = ^gi^i {cria + (7~o)) , (1) 

i 

where g is the coupling constant, jii is the projection of the 
transition dipole moment of a single chromophore to the po- 
larization of the cavity mode, erf are rising and lower opera- 
tors for chromophore excitations, and o) , a denote the cav- 
ity photon creation/annihilation operators. See the supple- 
mentary information (SI) for additional details concerning the 
model. The rotating wave approximation utilized in Eq. ([T]) is 
valid for the coupling regime we consider here. For the sake 
of simplicity we analyze photon statistics in a single polariza- 
tion mode in the cavity. In general, there are two degenerate 
polarization modes. However, the interaction between them 
is mediated by the LHC and the associated effects in the pho- 
ton statistics are of higher order in the LHC-cavity coupling as 
compared to the single photon emission by the LHC. The cou- 
pling coefficient g is given by ^ = a/SF' where = 1, V is 
the mode volume, and e is the permittivity of the medium be- 
tween the cavity mirrors. In the following, we assume £ = 2, 
which is approximately the dielectric constant of a solvent at 
optical frequencies. The cavity properties can be described 
by the quality factor Q and the mode volume V. While the 
former parameter determines the cavity photon escape rate 
as n = uJc/Q, the latter one controls the coupling with the 
LHC. Experimentally achievable values of the quality factor 
can be as large as Q = 10^ - 10^ l^il^El- The funda- 
mental limit of the mode volume is determined by the photon 
wavelength Vf = (A/2)^ and is dependent on the dielectric 
constant of the medium. While the mode volume in photonic 
crystal cavities is about to the best of our knowledge 

the state-of-the-art mode volume, demonstrated for a Fabry- 
Perot cavity design 1331 , which can be used for molecular 
spectroscopy, is V^xp ^ 5.5A^. This value corresponds to 
g « 0.015 cm~^/D coupling between the LHC exciton and 
the cavity vacuum state for Avacuum = 800 nm. A more op- 
timistic estimation for the coupling constant obtained for the 
fundamental limit of the mode volume, i.e. Vf ^ 0.023 /im^, 
is ^th ^0.1 cm~^/I^. In the following, we analyze the re- 
sults corresponding to the latter case in order to enhance the 



fundamental issues of molecular spectroscopy in a cavity we 
want to stress, while we also discuss (in SI) the state-of-the-art 
case to answer the question of which effects can be observed 
at the present stage of technological advances. As an exam- 
ple of an LHC, we choose the Fenna-Matthews-Olson (FMO) 
pigment-protein complex involved in the natural photosynthe- 
sis of green- sulphur bacteria. The FMO subunit has seven lf34ll 
strongly coupled bacteriochlorophyll (BChl) molecules (sites) 
d O . The Hamiltonian of the complex can be written as 

where only lowest electronic transitions of the molecules are 
taken into account, ujj is the site energy, and Vj^i denotes the 
coherent coupling between the corresponding sites. The sin- 
gle exciton transitions in the FMO complex are estimated to 
be in the range of 12150 - 12750 cm"^ |35, 36 1. The exci- 
tation dynamics in the FMO-cavity system can be described 
by a general Liouville-von Neumann equation for the density 
matrix 

p = -i[H,p\^C{p), (3) 

where the total Hamiltonian H includes the FMO Hamiltonian 
Eq. ([2]), the cavity Hamiltonian i^ph = uJcO^a, the interactions 
between FMO and the cavity Eq. ([T]), and also the coupling 
between the FMO and an external field E{uJi^t) given by 

Hi{t) = -^fiiE{uji,t)at ^h.c. (4) 

i 

The dephasing and relaxation channels in the isolated FMO 
system impose an irreversible dynamics of an initially cre- 
ated exciton state. The simplest way to introduce these noise 
processes into the dynamics is using a Lindblad form for 
the term C{p) where the environmental noise is described by 
the site energy fluctuations (phenomenological Haken-Strobl 
noise) |37| with dephasing noise strength 7. This analysis 
will bring out the basic principles of our proposal. Exten- 
sions to non-Markovian effects which can arise from strong 
coupling and/or the form of the environment spectral function 
ifTTl [T2I [38H411 as well as other cavity setups involving addi- 
tional laser fields, the role of multiple excitations, and other 
generalizations will be investigated in a forthcoming paper. 

Results and discussion- In Fig. [2j we show the mean 
photon number in the cavity, in the stationary state, as a 
function of the driving field and the cavity mode frequencies 
together with the electron excitation spectrum of FMO. This 
2D map is obtained 45 ps after a strong laser driving field 
of intensity / ^ 110 KW/cm^ has been turned on. The 
estimated values of electronic transition dipoles of the FMO 
chromophores are in the range of 3 — 14 D (see SI for the 
details of the model). Thus, the maximal coupling energy 
between the driving field and the excitons is of the order of 
1 — 7 cm~^ = 30 — 210 GHz and the exciton-cavity coupling 
is about 5 times smaller. The cavity photon population 
saturates to a stationary regime on 20 ps timescale showing 
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FIG. 2: Stationary mean photon number of the cavity mode as a 
function of the cavity (Q = 10^) resonance frequency ujc (vertical 
axis) and the pump laser frequency ui (horizontal), both in units of 
cm~^, but shifted by 12195 cm~^ ^ 820 nm, in the case of de- 
phasing rate 7 = 10 cm~^. The diagonal peaks correspond to the 
FMO exciton energies, while the off-diagonal peaks are related to 
energy transfer between different excitons. The asymmetry of the 
2D spectrum, with respect to the main diagonal (uji = cjc), is due to 
the different strengths of the laser and the cavity fields. Top Inset: 
computed spectrum of electronic excitations the FMO complex. The 
peaks corresponding to the FMO exciton energies are broadened by 
the 7 = 10 cm~^ linewidth. The red box is analyzed in Fig. [s] 



a set of peaks that are clearly associated with the exciton 
frequencies of the FMO complex. Indeed, in analogy to 2D 
spectroscopy lflj4ll36L diagonal peaks are in correspondence 
of the FMO exciton energies, and off-diagonal features 
appear due to energy transfer between different exciton states. 
The diagonal line in the spectra is due to both a coherent 
and an incoherent transfer of the photons from the laser field 
into the cavity. The former process is similar to the Raman 
scattering and goes through coherence pathways without 
generating an exciton population in LHC, while the latter 
process can be considered as a light absorption by the LHC 
with the following spontaneous emission of a same-frequency 
photon into the cavity. The horizontal lines in the spectra 
(off-diagonal peaks) are due to population transfer between 
different excitons. Weaker optical fields would result in 
a similar set of peaks, but the stationary regime would be 
achieved on a longer timescale. For the particular case shown 
in Fig. [2} the noise strength is 7 = 10 cm~^ ~ 2 ps~^ 1311421. 
Qualitatively, this picture remains valid for the values of 7 
within the range 1 — 100 cm~^ - see SI for more data. In 
all our simulations, we include an isotropic averaging of the 
measured quantities over completely random orientations 
of FMOs in the cavity. In presence of an oriented sample 
these quantum features would be further enhanced. Stronger, 
albeit (biological)-structure preserving, laser driving fields. 
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FIG. 3: Mean photon number of the stationary cavity mode as a func- 
tion of the cavity (Q = 10^) resonance frequency cjc (bottom axis) 
and the pump laser frequency ui (right axis), both in units of cm~^ 
- but shifted by 12195 cm~^ 820 nm - in the finely tuned range 
[245, 265] cm~^ (red box in Fig. [2]). One can clearly observe an 
anti-crossing effect. 



can result in a fine structure of a diagonal peak. This effect is 
similar to a level anti-crossing or a hole burning - see Fig.|3] 
Indeed, the relatively strong coupling with the pump laser 
field leads to a splitting of the exciton level into two dressed 
states corresponding to two new peaks in the mean photon 
number 2D map. In order to characterize quantum properties 
of the generated cavity photon state in its stationary state, we 
compute the second-order photon coherence function fl3l 




0.2 0.4 0.6 0.8 1.0 

FIG. 4: Second-order coherence function at time delay zero g^'^^ (0) 
of the stationary cavity mode as a function of the laser field coupling 
(horizontal axis) and the cavity (Q = 10^) resonance frequency cjc 
(vertical) - in units cm~^ - for pump laser frequency uoi — —uJc + 
510cm~^. The black line corresponds to g^'^^ (0) = 1. 
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FIG. 5: Photon mode population (rescaled by a factor 10 and 
second-order coherence function at zero time delay g^"^"^ (0) of the 
stationary cavity mode as a function of the dephasing rate (cm~^), 
for a diagonal (circles) and off-diagonal (square) peak in the mean 
photon number 2D map in Fig. |2] 
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at zero time delay g^'^\0) = ^j^^ - see Fig. 
general, we obtain a non-classical photon state geiieration 
(^•^^^(0) < 1) for most of the laser and cavity frequencies 
except of the diagonal line with Uc = uji where a coherent 
state {g^'^\^) = 1) and a thermal state (^"^^^(O) > 1) can be 
observed. More data on ^'^^^(0) in different regimes can be 
found in SI. Moreover, it turns out, as in Fig. [5] that this 
quantity ^*^^^(0) depends on the amount of dephasing noise 
in the dynamics and this might become a tool to estimate 
the strength of interaction between the biological molecule 
and the external environment. Notice also that the photon 
mode population is more sensitive to the dephasing noise, 
as compared to ^*^^^(0). In addition, the new scheme we are 
proposing allows us also to generated non-classical states of 
light by means of the coherent exciton dynamics. Finally, as 
possible signature of polariton formation, we quantify the 
amount of quantum correlations between the FMO sample 
and the cavity mode, as measured by the logarithmic nega- 
tivity [44 1 . It turns out that this quantity reaches a maximum 
at about 800 fs and around the main two diagonal peaks in 
the 2D map in Fig. |2]- see Fig. [6] Note that the lifetime of 
these correlations is around 2 ps~^, which is the time scale 
of the fastest decay process due to the presence of dephasing 
noise. Therefore, the coupling between the quantum cavity 
and the confined sample leads to the creation of non-classical 
correlations between these two systems. A different scheme 
towards entangling two FMO samples (entangled polaritons) 
in separate cavities is also discussed in SI. 

Conclusions and outlook- We have extended the rich 
physics of cavity quantum electrodynamics to biological 
molecules, in particular light-harvesting complexes involved 
in natural photosynthesis. As an example, we have consid- 
ered a sample of the Fenna-Matthews- Olson (FMO) pigment- 
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FIG. 6: Measure of quantum correlation (log-negativity) between the 
cavity mode and the FMO sample after ^ 800fs, as a function of 
the cavity (Q = 10^) resonance frequency uoc and the pump laser 
III frequency oji, both in units of cm~^, but shifted by 12195 cm~^, in 

Inset: Time evolution of 



the case of dephasing rate 7 = 10 cm 
the log-negativity measure for ujc = (^i 
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protein complex inside an optical cavity and driven transver- 
sally by a laser field. Our main result is that the emission 
spectrum from the coupled FMO-cavity system reflects coher- 
ent energy transfer into the cavity through delocalized exciton 
states, as well as exciton dynamics within the FMO complex. 
Moreover, we have found that a strong laser field driving the 
exciton dynamics in the protein complex can bum a hole in 
the emission spectrum resulting in an additional structure of 
resonance peaks. Finally, the generation of quantum states of 
light inside the cavity, due to the interaction with the confined 
biological complexes, was also observed. We believe that this 
novel spectroscopic tool can efficiently probe both static and 
dynamical (system-environment) properties of natural and ar- 
tificial photosynthetic structures and provide additional exper- 
imental support for the coherent dynamics unveiled by non- 
linear spectroscopy experiments. This combined evidence is 
expected to shed further light towards the elucidation of the 
role that quantum coherence may play in the remarkably ro- 
bust energy transport phenomena involved in natural photo- 
synthesis. 
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C02-01), the Basque government (project IT472-10), the 
SFB/TRR21, the EU Marie-Curie Programme, and the 
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THE MODEL 

In the site basis, the Hamiltonian of the FMO pigment- 
protein complex, Hs (in the main text), is given by 
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where the diagonal elements are the site energies - shifted 
from the base line c^base = 12195 cm~^ corresponding to a 
wavelength of = 820 nm - while the off-diagonal elements 
are the inter-site coupling rates (all numbers are given in units 
of cm-i = 1.988865 • 10"^^ Nm = 1.2414 lO""^ eV). 
The off-diagonal terms of the Hamiltonian were calculated 
in a dipole-dipole approximation and the site energies were 
taken from the Poisson-Boltzmann quantum chemistry model 
of Ref. Q. The intensities of the electronic transitions in the 
computed absorption spectrum (top inset in Fig. 1 of the main 
text) are proportional to l/i^p, and the spectrum is averaged 
over different spatial orientations of the FMO complex. The 
dissipation and dephasing caused by the surrounding environ- 
ment can be described by the following local Lindblad terms: 



^r,/2[-{a+a-,p} + 2a-pa+] , (5) 



with Tj and being the dissipative and dephasing rates at 
the site j, respectively, (a~) being the raising (lowering) 
operators for site j. In the following, we choose for simplicity 
uniform dephasing rates, i.e. equal and so simply labeled 
as 7. Moreover, here we neglect the presence of dissipation 
inside the FMO complex, because the 1 ns excitation lifetime 
f2\ is much longer than the time scale we look at. Let us point 
out that some dissipation is also present indirectly through the 
interaction with the cavity which has a damping channel, de- 
scribed by the Lindbladian term: 



^d^ssip) = rc/2 (- {ata, p} + 2apat) 



(7) 



where Tc being the damping rate of the cavity modes due 
to some dissipative processes, i.e. Fc = ^ with Q being 
the quality factor of the cavity. For a self-consistent model 
we use the same transition dipoles to construct interactions 
between BChl molecules in the FMO Hamiltonian and also for 



the coupling with the cavity and external fields. The positions 
and the relative directions of the dipoles are extracted from the 
structure of FMO complex | 3 1. The positions Ri of the 7 sites 
(BChl), taken as a middle point between 4 Nitrogen atoms in 
each BChl, are 



( 53.08 
56.04 
49.57 
38.81 
34.15 
41.44 



58.26 
54.79 
44.77 
42.25 
47.78 
47.82 



20.64 \ 

32.40 

45.39 

43.06 

31.26 

22.61 



(8) 



\ 47.53 43.95 33.22 / 

The directions the three components (in some reference 
frame) of the 7 induced-transition dipoles pi are 



[Pi] 
Po 



\ 



-0.026 

-0.752 

-0.935 

-0.001 

-0.739 

0.859 

0.176 



0.286 
0.601 
0.061 
0.393 
0.672 
0.371 
-0.042 



(9) 



/ -0.026 0.286 -0.958 \ 
-0.271 
0.349 
-0.919 
0.048 
-0.353 
-0.983 J 

where the absolute values of the transition dipoles are taken 
as a phenomenological parameter /io = 6 D, which accounts 
for effects of screening and induced charges |T| . The exciton 
energies corresponding to the Hamiltonian ([5j are £^i_7 = 
[-31.0, 129.4, 148.7, 254.6, 310.1, 394.5, 518.6] cm-^ 
and the absolute values of the electronic transition 
dipoles of the FMO chromophores are l/ii-rl = 
[6.4,13.8,3.7,11.3,7.2,4.0,5.6] D. 

Finally, we assume a continues-wave laser pumping exciton 
transitions in the LHC. The laser field is coupled to the cavity 
modes through the trapped bio-sample only. The value of the 
external (pump) field can be estimated as 



E = 



(10) 



where / is the laser field intensity, c is the speed of light, and 
So is the vacuum permittivity. For example, the intensity of 
a CW laser field used in Raman spectroscopy is about 5 mW. 
In these experiments the laser beam is focused on a spot of a 
size of 5 /im^. This corresponds to / = 100 kW/cm^ or the 
field E = 27 ■ 10^ V/m = 0.45 cm-7^. This value is much 
smaller than the peak field value used in femtosecond pulse 
lasers. However, a CW would result in a heat accumulation, 
thus stronger fields could damage a sample. 

EXTRACTING ADDITIONAL INFORMATION 

In this section we investigate other important features of the 
LHC dynamics that we can extract measuring statistical prop- 
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FIG. 7: Stationary mean photon number of the cavity mode as a 
function of the cavity (Q = 10^) resonance frequency ujc (vertical 
axis) and the pump laser frequency ui (horizontal), both in units of 
cm~^, but shifted by 12195 cm~^ ^ 820 nm, for dephasing rates 
7 = 1, 50, 100 cm~^ (from top to bottom). 




FIG. 8: Average photon number of the stationary cavity mode as a 
function of the laser field coupling (horizontal axis) and the cavity 
(Q = 10^) resonance frequency uJc (vertical) - in units cm~^ - for 
pump laser frequency uoi — —uJc + 510cm~^, which corresponds 
to an cross-diagonal cut of a diagonal peak (red box in Fig. 2 of the 
main manuscript). 




FIG. 9: Stationary mean photon number of the cavity mode as a 
function of the cavity (Q = 10^) resonance frequency cjc (vertical 
axis) and the pump laser frequency uji (horizontal), both in units of 
cm~^ - shifted by 12195 cm~^ - for dephasing rate 7 = 10 cvnT^ . 



the present stage of technological advances. The behaviour of 
the anti-crossing effect withing a diagonal peak as a function 
of the laser field coupling can be observed in Fig. [8] 



Inhomogeneous broadening: static disorder 



erties of the outcoming light (probe) mode. The stationary 
mean photon number of the outcoming light mode is shown 
for different values of dephasing rate 7 (in Fig. [7]), also for a 
lower-quality {Q — 10^) cavity (in Fig. [9]), and for a weaker 
(state-of-the-art) cavity coupling rate (in Fig. [T0| according to 



Here, we show the effect of static disorder, i.e. disorder in 
the site energies, on the mean photon number 2D spectra, av- 
eraged over full random orientation disorder as well (as done 



in the rest of the manuscript) - see Fig. 1 1 for a disorder of 50 
cm~^. As observed in 2D spectroscopy, it lead to the presence 
of inhomogeneous broadening for the diagonal peaks. 
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FIG. 10: Stationary mean photon number of the cavity mode as a 
function of the cavity (Q = 10^) resonance frequency Uc (vertical 
axis) and the pump laser frequency uji (horizontal), both in units of 
cm~^ - shifted by 12195 cm~^ - in the case of dephasing rate 7 = 
50 cm~^, weaker pump field E = 0.01 cm~^/D, and weaker cavity 
coupling rate g — 0.01 cm~^ / D. 
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FIG. 1 1 : Stationary mean photon number of the cavity mode as a 
function of the cavity (Q = 10^) resonance frequency ujc (vertical 
axis) and the pump laser frequency uoi (horizontal), both in units of 
cm~^ - shifted by 12195 cm~^ - in presence of an inhomogeneous 
broadening, induced by a static disorder in the site energies, set to 50 
cm~^, in the case of dephasing rate 7 = 10 cm~^ . 



Second-order correlation function 

Here, we analyze the time evolution of the second-order co- 
herence function at time delay zero ^'^^^(0) of the stationary 
cavity mode (defined in the main text) for diagonal and off- 
diagonal peaks - see Fig. 



12 Although the time scale, we 



consider here, is not feasible to be investigated with the cur- 
rent photon detector technology, however it could be in a near- 
future or even nowadays with other natural or artificial light- 
harvesting systems where this behavior takes place in a longer 




4 6 
time (ps) 

FIG. 12: Second-order coherence function at time delay zero g^'^^ (0) 
of the cavity (Q = 10^) mode versus time (ps) for different values of 
dephasing rates and for diagonal and off-diagonal peaks, in the case 
of dephasing rate 7 = 10 cm~^ . 



time due to different system parameters. Furthermore, in Figs. 
13 we show a full-range contour-plot for the second-order co- 
herence function at time delay zero ^'^^^(0) of the stationary 
cavity mode as a function of the cavity and pump frequen- 
cies, for different cavity parameters. We observe observe a 
formation of a quantum photon state {g^'^^ (0) < 1) almost ev- 
erywhere, except along the diagonal line where coherent and 
thermal states (^*^^^(0) > 1) can appear. 



TRAPPING LHCS INSIDE THE CAVITY 

In our analysis it is important to trap the biological sample 
inside the cavity. Although several techniques are know in this 
respect, in this section we show, by a simple and elementary 
argument, the feasibility of such process by using an external 
electric field. In particular, let us apply an additional constant 
electric field Eq, polarized along the axis e with carrier fre- 
quency uou interacting with the FMO complex as the pump 
laser field. The corresponding interaction energy is 



E 



I Mi 



^0 



UJi - Ldi 



(11) 



and shows a minimum by varying the orientation of the polar- 
ization axis e in terms of the two angles 6 and (j) with respect 
the site-1 dipole moment, i.e., for 6> ~ 1.75 and (j) ^ 2. In 
particular, we consider a spatial profile for the electric field 
given by a Gaussian function peaked in the center of the 



cavity as Eo{x) 
Fig 



Eoe- 



where A = 800 nm. In 



14 we show the trapping potential A(x) and the cor- 



responding trapping force as a function of the position x, in 
the case of i9 1.75, (j) 2, uji = -5000 cm"\ and 
^0 = 5 cm-^ - 30 • 10^ J/{C m). The trapping 
force is about 20 times larger than the gravitational force to 
which a single FMO is subjected, assuming that it has a mass 
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FIG. 13: Second-order coherence function at time delay zero g^^^ (0) 
of the stationary cavity {g — 0.1 cm~^/D and Q = 10^ (Top) 
and Q = 10^ (Middle); Bottom: g = 0.01 cm"V^' ^ = 
0.01 cm" "VA andQ = 10^) as a function of the cavity resonance 
frequency cjc (vertical axis) and the pump laser frequency uji (hori- 
zontal), both in units of cm~^ - shifted by 12195 cm~^ - for dephas- 
ing rate 7 = 10 cm~^ . The black lines refer to g'^^^ (0) = 1. 



of around 80 kDa ~ 15 10~^^Kg (including the protein scaf- 
folding) |4|. In this regime, the amount of electronic excita- 
tion in the FMO complex is very small - see Inset in Fig. 14 
Moreover, we have found that the trapping force can be even 
several order of magnitudes larger than the gravitational one, 
allowing a very tiny amount of excitation in the FMO system 
induced by the presence of the laser irradiation. Therefore, 
following this simple argument, it seems that it is possible to 



FIG. 14: Trapping force (in units of gravity force) vs. position (nm) 
for a single FMO complex, when subject to a constant laser field with 
a spatial Gaussian profile peaked in the center of the cavity, ^0 [x) — 
^^g-x2/(2A) ^-^j^ ^0 = 5 D-^ cm-^ - 30 • 10^ J/{C m), car- 
rier frequency set to ui = —5000 cm~^ and polarization orientation 
given by ^ ^ 1.75, ~ 2. Top inset: trapping potential (in cm~^) 
as a function of the position (in nm). Bottom Inset: in the presence 
of such laser irradiation (even applied for about 1 ps), the amount 
of total excitation in the FMO system is essentially vanishing, i.e. 
smaller than 3 • 10~^. 



trap a sample of FMO complexes inside the cavity, without 
perturbing the system, i.e. leaving it in its ground state, and 
then experimentally apply a driving laser field as described 
above. 



ENTANGLING BIO-SAMPLES IN SEPARATE CAVITIES 



Finally, we propose a possible scheme (in Fig. [T5]) for a 
probabilistic creation of quantum correlations, i.e. entangle- 
ment, between two FMO samples confined in spatially sep- 
arate cavities, following Ref. |5 1 where such protocol was 
introduced for trapped ions. In other terms, after applying the 
pump laser fields to the two separate samples, the outcoming 
photons from the cavities are mixed on a 50/50 beam-splitter 
and one applies a projection of the two-photon state into a Bell 
state of the form l/>/2(| click, no click) + | no click, click)). 
Then, after the projection, the composite system of the two 
LHC samples has a chance to be left in an entangled state 
with some finite probability. Therefore, we measure the en- 
tanglement (by log-negativity) between the two samples, after 
the projection for the two-cavity mode into the Bell state. This 



quantity is shown in Fig. 16 We find that one can probabilis- 
tically entangle the two biological samples (entangled polari- 
tons) with a non- vanishing probability. This scheme may, for 
instance, allow us to extract more information on the quan- 
tum performance of these biological molecules, and perhaps 
pave the way for new tests about the foundation of quantum 
mechanics in terms of biological systems where quantum fea- 
tures seem to play an important role. 
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FIG. 15: We consider a setup of two equal spatially separate optical 
cavities in which two LHC samples are confined, while subjected to 
two similar pump laser fields. Two auxiliary photodectors are also 
introduced to detect photons, leaking out of the cavities, after being 
mixed on a 50/50 beam-splitter. 




FIG. 16: Amount of quantum correlations (measured by logarith- 
mic negativity) between the two FMO samples trapped in two spa- 
tially separate cavities, after ^ 140fs, in the case of dephasing rate 
7 = 10 cm~^. Inset: Entanglement vs. time (ps) for Uc — (jOi — 
200 cm"\ 
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